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Langmuir  Shock  Waves 


Martin  V.  Goldman*  and  Dwight  R.  Nicholson** 
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Boulder,  Colorado  80309 


ABSTRACT 


Intense  Langmuir  waves  in  a uniform,  one-dimensional 
plasma  are  described  by  a nonlinear  Schrodinger  equation. 
When  the  effect  of  linear  Landau  damping  is  included,  we 
find  steady  state  shock  wave  solutions.  For  an  infinite 
plasma,  frequency  determines  shock  amplitude  and  (subsonic) 
speed.  For  a semi-infinite  plasma,  the  shock  is  stationary; 
its  frequency  determines  its  amplitude,  and  the  energy  flux 
through  the  boundary  determines  the  spatial  extent  of  the 
field. 
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INTRODUCTION 


It  is  well  known  that,  under  certain  conditions, 
intense  Langmuir  waves  characterized  by  electric  field 
[e  (x,  t)  exp  (-ia)p0t)  + c.c.]  and  the  associated  ion  density 
perturbation  n^(x,t)  are  described  in  time  and  in  one 
spatial  dimension  by  the  dimensionless  equations 


i - n.E  + iL(E)  = 0 

9t  9x2  1 


, ra.  3 n. 
2 _i  x 

4 m . 2 
e 3t 


92|fi|2 


where  we  have  included  a linear  damping  operator  L (E)  in 

(1)  and  m.  and  m are  the  ion  mass  and  electron  mass  re- 
l e 

spectively.  Physical  variables  are  obtained  from  the 

dimensionless  variables  in  (1)  and  (2)  as  follows: 

1/2  1/2 

distance  is  3 Ax,  where  A = (T  /m  ) is  the  electron 

e e e e 

Debye  length;  time  is  equal  to  2t/u>p0,  where  u)po  is  the 

unperturbed  background  electron  plasma  frequency;  the 

. - 1/2 
electric  field  is  equal  to  E (x,  t)  (4n»eTe)  co^/e,  where 

e is  the  absolute  value  of  the  electron  charge;  ion 

density  perturbation  is  n^nQ,  where  nQ  is  the  unperturbed 

ion  density;  the  ion  temperature  T^  = 0 and  E(x,t)  and  n^(x,t) 

vary  on  a time  scale  long  compared  to  u 


An  extensive  bibliography  of  work  on  these  equations  can 
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be  found  in  Nicholson  and  Goldman,*  who  study  the  temporal 
damping  of  soliton-like  initial  conditions.  Appendix  A 
of  the  present  paper  provides  a list  of  the  assumptions 
which  enter  the  derivation  of  (1)  and  (2). 

Suppose  there  is  a frame  (where  the  ion  disturbance 
is  stationary)  moving  with  velocity  V with  respect  to  the 
laboratory.  If  z = x-Vt,  we  can  solve  (2)  for  n^;  defining 
E(x,t)  = E(x,t)[l  - V2(3nr/4me)]  *^2,  Eq.  (1)  becomes 


M - iV 
at 


3E 

az 


a2E 

az2 


+ |ETE  + iL  (E)  = 


(3) 


1/2 

This  result  is  valid  for  speeds  |v|  <<  (4me/3nu)  ' or  in 

1/2 

physical  units,  speed  less  than  sound  speed  c = (T  /m. ) 

S 0 1 

Using  Eq.  (3) , Nicholson  and  Goldman*  studied  the 

damping  of  a solitary  wave  initial  condition.  This  work 

2 3 

is  analogous  to  the  work  of  Ott  and  Sudan  ' on  the  Korteweg- 
de Vries  equation. 

In  the  present  paper,  we  discuss  solutions  of  (3) 
which  exhibit  shock-like  behavior;  that  is,  the  field  at 
x -*  -<*>,  |E(x  = -°°,  t)  | has  a different  value  than  it  does 
at  some  other  point  in  space,  either  x = +®  for  the  case  of 
infinite  plasma,  or  x = 0 for  the  case  of  a semi-infinite 
plasma. 

Section  II  presents  an  analytic  approach  to  Eq.  (3). 


1 


In  Section  III  we  discuss  numerical  solutions  to  (3) 

with  an  infinite  medium  and  L(E)  representing  linear  Landau 
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damping;  here  we  find  that  the  shock  amplitude  implies 
shock  speed  and  vice  versa.  Section  IV  discusses  steady 
state  shock  solutions  with  zero  speed  for  a semi-infinite 
medium,  with  a constant  energy  flux  injected  at  the  boundary. 
Conclusions  and  related  work  are  discussed  in  Section  V. 

II.  ANALYTIC  CONSIDERATIONS 

We  look  for  a solution  of  (3)  of  the  form 

EU.t)  - p1/2(2)  ei°<2l+it!t  , ,4) 

1/2 

where  the  amplitude  p (z)  is  independent  of  time. 

Inserting  (4)  into  (3)  and  separating  real  and  imaginary 
parts  we  have 


1/2 

(In  Appendix  B we  discuss  Eq.  (5)  rewritten  in  terms  of  p 
instead  of  p;  this  is  more  instructive  for  some  purposes.) 
Equation  (5)  is  analogous  to  a particle  energy  equation, 
of  the  form 

\ p'2  + $(p)  = W (p,  z)  > 0 , (9) 

No  Damping 

2 

where  (1/2)  p*  is  the  "kinetic  energy",  <t>(p)  is  the 
pseudopotential,  and  W(p,z)  is  the  "total  energy";  and 
where  the  amplitude  p(z)  plays  the  role  of  "particle 
velocity"  and  the  position  z plays  the  role  of  "time." 

The  "total  energy"  W(p,z)  goes  to  zero  as  z -*■  -°°,  and  is 
identically  zero  in  the  absence  of  damping.  The  pseudo- 
potential 4>(p)  = p^  + 2p2[(v2/4)  - ft]  . 

For  zero  damping,  we  obtain  the  soliton  solution 
as  shown  in  Fig.  (la).  A "particle"  in  the  potential 
well  takes  an  infinite  "time"  z to  roll  from  the  initial 
"position"  p = 0 to  the  reflection  "position,"  p = p„,„, 

MAX 

and  an  infinite  amount  of  "time"  z to  roll  back  again. 


In  this  undamped  case,  Eq.  (6)  is  trivially  solved  to 
obtain  o(z)  = Vz/2 , and  the  pseudopotential  equation  (5) 


6 


can  be  analytically  solved  for  p(z).  After  transforming 
back  to  the  laboratory  frame  x = z+Vt,  we  have  the  well- 
known  class  of  single  soliton  solutions 

E(x,t)  = [2Q  - 


J\l/2 


t)  sech 


17  v2\1/2  1 

~ x/  (x'vt) 


or,  in 


exp^iQt  - 2 

terms  of  A = [ft  - (V2/4)] 


(9) 


1/2 


E (x,  t)  = 2^2A  sech  [a  (x-Vt)  ] exp^iA2t 


iY^t  iv 

4 2 


(10) 


) 


which  exist  for  0 < A2  <<  [l  - V2 (3/4) (m./m  )]  and 
| V | < (4me/3mi)1/2. 

With  the  inclusion  of  Landau  damping  in  (5) -(8),  it 
can  be  verified  numerically,  as  in  the  following  two  sections, 
that  the  total  energy  W(p,z)  decreases  slowly  (not  always 
monotonically)  with  "time"  z as  the  "particle"  rolls  in  the 
well,  yield  a shock  solution  as  shown  in  Fig.  (lb) . In 
Section  III  we  numerically  obtain  an  infinitely  extended 
shock;  the  shock  speed  is  determined  by  the  shock  amplitude. 

In  Section  IV  we  find  semi-inf initely  extended  shocks,  which 
are  stationary  and  have  an  energy  flux  through  the  boundary. 


FIGURE  lb 


Ill . 


INFINITE  SHOCK 


We  obtain  an  infinitely  extended  Langmuir  shock 
wave  by  treating  Eq.  (3)  numerically  with  the  inclusion  of 
the  linear  Landau  damping  operator^" 


- @r  / * 


elkz  E (k, t) 


exp  [ (-3/2) k2] 

W1 


dz ' (z-z ' ) E (z ' , t) 


oo 

/ dk  exp[(-3/2)k2Jsin  [k (z-z1 ) ] 


We  again  look  for  a solution  of  the  form 


E (z, t)  = p1/2(z)eia(z)+ifit  , 


(ID 


(4) 


with  p(z)  independent  of  time  (in  this  moving  frame). 
Numerically,  we  guess  a function  E(z,t)  to  insert  into  the 
damping  operator  L in  (11).  Specifically,  we  guess  the 
soliton  in  Eq.  (9) ; we  then  numerically  integrate  the 
(now  time- independent)  Eq.  (3)  from  z -*■  -°°,  taking  the 
soliton  form  of  the  solution  from  (9)  as  the  boundary 
condition  at  large  negative  z.  The  resulting  E(z,t)  is 


*-■— 


-a 
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then  inserted  into  the  damping  operator  (11),  and  the 
entire  procedure  is  iterated  until  convergence  is  obtained. 
(The  convergence  is  usually  quite  rapid. ) 

The  only  allowable  solutions  are  those  for  which 
a(z)  ♦ constant  as  z -*•  +°°,  and  p(z)  -*•  constant  as  z -*•  +°°. 

If  this  were  not  the  case,  then  the  Landau  damping  operator 
would  be  finite  over  an  infinite  region  of  space,  yielding 
an  infinite  amount  of  energy  dissipation  per  unit  time  and 
no  possibility  of  a steady  state  shock  solution. 

The  constancy  of  o(z)  and  p(z)  at  large  z implies 
that  there  can  be  no  energy  flux  from  large  z into  the 
region  of  Landau  damping.  To  establish  this  fact,  we  use 
energy  conservation  laws:  multiply  (3)  by  E*,  multiply 

the  complex  conjugate  of  (3)  by  E,  subtract  and  integrate 
over  all  space;  we  find  in  the  shock  frame,  using  the  fact 
that  E ( z -»  t)  -*■  0, 


-V|E| 


+ 2Im[E*OE/3z)]  I 

Z-*-+°o  [ Z-*-+°° 


oo 

/ 


= -2Re  I dz ' E*L (E)  . 


(12) 


The  second  term,  assuming  p(z)  constant  as  Z +oo  9 is 


2p(z)  (3o/3z) 


(13) 


Z+oo 


■ ! 


This  term,  representing  energy  flux  from  z -*•  +o°,  vanishes 


since  3o/3z|z^m  + 0 by  the  argument  above.  The  remaining 
two  terms  show  that  the  energy  dissipated  by  Landau  damping 
(right-hand  side)  is  exactly  balanced  by  — V | E | |z_^+0O»  which 

is  just  the  rate  at  which  electric  field  energy  is  disappear- 
ing from  a fixed  spatial  volume  of  the  laboratory  frame 
because  the  shock  is  moving  to  the  right  with  speed  V. 

Now  we  already  know  the  value  of  I E I I _ . = p in 

terms  of  the  pseudopotential  of  Fig.  (lb);  this  is  just  the 
point  where  d<J>(p)/dp  + 0 or  p#  = (4/3) [ft  -(V2/4)]  = <2/3)pmax- 
Thus,  (12)  becomes  an  implicit  expression  for  V; 

V = — Re 

2[fi-(V74)] 

in  which  it  must  be  remembered  that  E(z,t)  obtained  numeri- 
cally depends  on  the  value  of  V used  in  (3).  Equation  (14) 
can  also  be  obtained  from  (6)  and  (7)  by  setting 
o'(z-m-°°)  = o.  In  practice,  we  must  numerically  repeat  this 
calculation  with  different  values  of  V until  we  find  the 
one  which  makes  da/dz  | 2_*+00  ■*  0. 

The  results  of  this  procedure  are  shown  in  Fig.  2 
for  ft  = 0.3;  the  value  V = 0.35  is  found  numerically  to 
make  do/dzl  , ■*  0. 

The  energy  conservation  law  (12)  indicates  (it  may 
be  difficult  to  prove  this  rigorously)  that  for  each  shock 
amplitude  there  is  a shock  speed,  and  vice  versa.  For  a 


-5  O 

Distance  2 
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msE  cr '(Hr) 


-5  O )5 
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FIGURE  2 
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given  electron-ion  mass  ratio,  we  require  I V | <<  (4m  /3m  .)  ' . 

e l 

Thus,  there  will  be  a maximum  shock  amplitude  for  each  mass 
ratio,  above  which  the  speed  needed  to  satisfy  the  energy 
conservation  law  (12)  is  higher  than  allowed  by  the  theory. 

It  is  interesting  to  discuss  our  results  in  terms  of 

4 

the  result  of  Pfirsch  and  Sudan,  who  proved  that  electron 
Landau  damping  of  ion  acoustic  waves  does  not  lead  to  a 
shock  solution  of  the  Korteweg-deVries  equation.  Pfirsch 

4 

and  Sudan  showed  that  a necessary  condition  for  shock 

solutions  to  the  Korteweg-deVries  equation  is  that 

LIM[y(k)/k]  -*■  Q,  where  (y(k)}  are  the  eigenvalues  of  the 
k+0 

1 lry 

damping  operator  L(E)  with  eigenfunctions  {e  }.  For 
electron  Landau  damping  of  ion  acoustic  waves,  y(k)  <*  k and 
the  necessary  condition  is  not  satisfied.  Now  it  turns 

4 

out  that  the  development  of  Pfirsch  and  Sudan  carries  over 

exactly  for  the  case  of  the  present  paper;  here,  electron 

Landau  damping  of  Langmuir  waves  has  y(k)  « exp [(-3/2) k ] / | k | J 

which  indeed  satisfies  the  necessary  condition 

LIM[y(k)/k]  - 0. 
k-"0 

In  this  section  the  normalizations  entering  the  non- 
linear Schrodinger  equation  are  slightly  different  in  the 
high  amplitude  half-space  than  in  the  low  amplitude  half- 
space; e.g.,  there  must  be  a very  small  particle  flow 
velocity  in  the  lab  frame,  high-amplitude  half-space,  to 
supply  particles  to  the  region  where  the  shock  front  has 
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just  passed.  These  differences  in  normalization  are 
ignored,  consistent  with  the  derivation  of  the  basic 
equations  (Appendix  A)  and  with  previous  work  on  nonlinear 
Schrodinger  equation  shocks.5 


IV.  SEMI-INFINITE  SHOCK 

Let  us  now  consider  a semi-infinite  medium, 

-°°  < x < L.  We  again  look  for  shock  solutions,  but  here  we 
allow  an  energy  flux  through  the  boundary  at  x = L,  travel- 
ing to  the  left.  We  use  the  numerical  procedure  of  the 
previous  section,  modifying  the  spatial  integral  in  (11) 
so  that  the  limits  of  integration  are  z -*•  -00  and  z = L. 

The  results  are  qualitatively  similar  to  those  of  Fig.  (2) 
for  the  infinite  medium,  with  two  differences:  namely, 

that  the  right-side  of  the  graph  now  represents  the  physi- 
cal boundary  at  z = L,  with  boundary  conditions  at  that 
point  on  |E(z=L)|,  o(z=L),  and  their  derivatives;  and 
that,  in  contrast  to  the  previous  section,  there  is  no  re- 
quirement that  o'(z-»-L)  0,  since  this  derivative  is  indeed 

proportional  to  the  energy  flux. 

Numerically,  the  frequency  of  the  shock  determines 


1 


the  shock  amplitude,  and  the  length  of  the  system  determines 
the  energy  flux  required.  Physically,  we  would  like  to 
make  the  argument  in  the  opposite  direction.  That  is, 
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! 


given  a semi-infinite  system  with  a boundary  at  z = L,  we 
apply  an  electric  field  at  a given  frequency  and  with  the 
correct  amplitude  for  that  frequency;  and  we  adjust  the 
spatial  derivative  of  the  electric  field  so  that  there  is 
a certain  energy  flux.  Then  a response  moves  into  the 
plasma,  forming  a shock  wave;  when  it  has  penetrated  so 
far  that  the  energy  dissipation,  integrated  over  -<*>  < z 
< L,  exactly  balances  the  energy  flux  through  the  boundary, 
the  shock  stops  moving  and  we  have  a steady  state  which 
looks  qualitatively  like  that  of  Fig.  (2)  . 

If  the  above  scenario  actually  occurs,  it  may  have 
many  implications  for  laser  fusion,  microwave  heating  of 
magnetically  confined  plasma,  radio  wave  modification  of 
the  ionosphere,  and  many  other  physical  situations.  The 
question  of  whether  the  above  scenario  actually  occurs  is 
strongly  related  to  whether  the  shock  wave  is  stable.  We 
have  undertaken  the  stability  analysis,  which  will  be  dis- 
cussed in  a separate  publication.  The  evidence  for  sta- 
bility from  the  two  known  limiting  cases  is  ambiguous;  the 
high  amplitude  part  of  the  infinite  shock  of  the  previous 
section  is  known  to  be  modulationally  unstable,  while  the 
soliton  shape  itself  is  known  to  be  stable.  Thus,  a care- 
ful stability  analysis  is  required  for  the  semi-infinite 
shock. 


| 


•i 
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V.  CONCLUSIONS  AND  RELATED  WORK 


In  an  infinite  medium,  we  have  found  that  the 
combination  of  electron  Landau  damping  and  ponderomotive 
force  nonlinearity  leads  to  the  existence  of  shock- like 
solutions,  traveling  to  the  right  (toward  the  high  ampli- 
tude region)  in  order  that  the  dissipation  energy  can  be 
balanced  by  a disappearance  of  electric  field  energy  in 
the  system.  There  is  a one-to-one  correspondence  between 
shock  amplitude  and  shock  speed;  for  a given  electron-ion 
mass  ratio  there  is  a maximum  allowable  shock  amplitude. 

A semi-infinite  shock  is  found,  stationary  in  the 
lab  frame,  with  energy  flux  balanced  by  dissipation. 

Our  work  differs  from  that  of  Mima,  Nishikawa,  and 

Ikezi.^  We  treat  electron  nonlinearity  and  electron  Landau 

damping;  our  shock  speed  is  0 < |v|  < c in  physical  units. 

* s 

They  treat  electron  nonlinearity,  ion  nonlinearity  (v  • Vv) , 
and  reflected  ions;  their  shock  speed  is  limited  to 
| V | ~ cg  in  physical  units.  We  expect  that  very  close  to 
the  ion  sound  speed,  ion  nonlinearities  will  become  impor- 
tant in  our  work,  and  the  treatment  may  need  modification. 


APPENDIX  A:  ASSUMPTIONS  IN  DERIVING  THE  NONLINEAR 

SCHRODINGER  EQUATION 

The  basic  equations  (1)  and  (2)  are  derived  from 
the  force  equation,  continuity  equation,  and  Poisson's 
equation  under  the  following  assumptions: 

(i)  All  quantities  can  be  characterized  as  high 
frequency  (2.)  or  low  frequency  (s)  . 

(ii)  Quasineutrality:  (a)  The  oscillatory 

distance  of  an  electron  in  the  high  frequency  electric 
field,  xQSC,  must  be  much  less  than  the  low  frequency 
scale  length,  LS;  *qSC  <<  Ls.  (b)  We  require  Ae  <<  L, 
where  X is  the  electron  Debye  length. 

(iii)  Weak  turbulence:  |E  I <<  n T leads  to 

1 o e 

s & 

<<:  A and  to  15,1  =<$n.  <<  n ' where  E is  the  high  fre- 

OSC  e e l o J 

quency  electric  field,  nQ  is  the  background  density,  and 

<5n  ■ is  the  (electron,  ion)  density  disturbance. 

6 # 1 

. S 

(iv)  Where  Vg  is  the  low  frequency  speed,  is 

l 

the  ion  speed,  and  vg  is  the  high  frequency  electron 
speed,  veS  ^ . 
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APPENDIX  B:  ALTERNATE  PSEUDOPOTENTIAL 


In  Section  II  the  pseudopotential  equation  (5) 

uses  p ( z ) as  a variable.  If  we  rewrite  that  equation  in 

1/2 

terms  of  n(z)  = p (z)  , we  obtain 


+ J n4  + \ n2(^~  - «)  « w(n,z)  , 


ao 

w(n,z)  = i |dz' [n2(z')]  ’ ri2  (z  • ) - ^^-1,  (A2) 

J L n (2 • > J 


and  I and  R are  defined  in  (7)  and  (8) . We  have  found  that 
the  right-hand  side  of  the  former  pseudopotential  equation 
(5)  is  not  a monotonically  decreasing  function  of  z;  it  is 
interesting  to  speculate  whether  the  right-hand  side 
(W(n , z) ) of  the  new  pseudopotential  equation  (Al)  is  a 
monotonically  decreasing  function  of  z.  Numerical  work  to 
answer  this  question  is  in  progress. 
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